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We report on electrical properties of the two sensing proteins: bacteriorhodopsin and rat olfactory 
receptor OR-I7. As relevant transport parameters we consider the small-signal impedance spectrum 
and the static current-voltage characteristics. Calculations are compared with available experimental 
results and the model predictability is tested for future perspectives. 

PACS numbers: 87.14.ep, 87.15.hp, 87.15.fh 



Sensing proteins are of relevant interest because of the 
fundamental role they play in living environments at a 
cellular level Recently, much work has been devoted 
to the investigation of the G protein-coupled receptors 
(GPCRs), a large family of proteins sensitive to the cap- 
ture of single or a few specific molecules (ligands) and a 
similar protein, the bacteriorhodopsin (bR), a light sen- 
sitive protein whose action is related to a proton pump. 
The activation of sensing proteins is a quite complex 
mechanism which in any case starts with a conforma- 
tional change. The capture mechanism and also their 
3D (tertiary) structure is quite the same for all of them, 
despite of the very different nature of the possible lig- 
ands (from photons to neurotoxins, from smells to hor- 
mones). The attention on these receptors is very high 
and the attempts of characterization cover different pro- 
tein features. On this respect, we mention recent exper- 
iments on the current-voltage (I-V) characterization of 
bR nanolayers when sandwiched between metallic con- 
tacts 0, Q , and when using an atomic force microscope 
technique [J]. On the same subject, recent experiments 
[1] were performed on the rat 0R-I7 anchored on a func- 
tionalized gold substrate with the technique of molecu- 
lar self-assembly. These measurements 0,1^ showed the 
possibility to monitor the protein sensing action, and the 
related conformational change, by means of the modifi- 
cation of the electrochemical impedance spectrum (EIS) 
in the presence of a controlled flux of specific odorants 
(heptanal, octanal). 

To our knowledge, the above experimental results did 
not frame into a unified theoretical scheme, neither from 
a macroscopicpoint of view Q, nor from a microscopic 
point of view Si. As an alternative to these approaches 
we developed a coarse grained model which connects elec- 
trical and morphological protein modifications by using 
an impedance network protein analogue (IPNA) [10l - [l3| . 
Here, the model is further validated on available exper- 
iments concerning the small signal impedance spectrum 
of rat 0R-I7 and its predictability is tested by a compar- 
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FIG. 1: (Color online) Contact map of the native and acti- 
vated state of bacteriorhodopsin (fig. la) and OR-I7 (fig. 2b). 
Fig. la: the native state (black full circles) vs. the activated 
state (red full triangles) for ii = 6 A . Fig. lb: the native state 
(violet full squares) vs. the activated state (magenta full tri- 
angles) for ii = 12 A . 



ative investigation on bR and rat 0R-I7. 

To set up the IPNA, the protein tertiary structure, in 
the native and activated state, have to be known. This 
information is available in the protein data banks, al- 
though for few proteins [13] ■ When these data have been 
gained, the protein in a given state can be mapped into 
a topological network (graph) [l^ [l^ , each node of the 
graph corresponding to a Cq,, taken as the single inter- 
acting center. When two Cq,s are closer than an assigned 
cut-off distance, R, a link is drawn between the corre- 
sponding nodes [l^]- Thus, by fixing the value of R, 
we fix the complexity of the graph. For bR we take the 
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FIG. 2: (Color online) Fig2a: Nyquist plot of bR for the PDB 
entries 2NTU, the native state, (continuous line) and 2NTW, 
the activated state, (dashed line). Plots are calculated for R — 
6 A (black) and 12 A (magenta), and values are normalized to 
the zero-frequency value of the impedance pertaining to the 
native state. Fig2b: Normalized real (blue) and imaginary 
(orange) parts of the total impedance of bR versus the angular 
frequency. The data are calculated for an interaction radius 
R = gA . Continuous lines refer to the native state, dashed 
lines to the activated state. 
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FIG. 3: (Color online) Fig2a: Nyquist plot of rat OR-I7, in 
the native state, (continuous line) and the activated state, 
(dashed line) . Plots are calculated for R = Q A (black) and 12 
A (magenta) , and values are normalized to the zero-frequency 
value of the impedance pertaining to the native state. Fig2b: 
Normalized real (blue) and imaginary (orange) parts of the 
total impedance of rat OR-I7 versus the angular frequency. 
The data are calculated for an interaction radius R= 12 A . 
Continuous lines refer to the native state, dashed lines to the 
activated state. 



couple of PDB entries 2NTU (native) and 2NTW (acti- 
vated) Q- For rat 0R-I7 we used the GPCR bovine 
rhodopsin as a template and the MODELLER software 
to construct, by homology, the tertiary structure for both 
the native and activated states [l2| . 

A comparative analysis between the native and acti- 
vated states of bR and rat OR 17, is given by their con- 
tact maps as reported in Figs, la and lb, respec- 
tively. For a fixed cut-off radius, these maps provide a 
reduced representation of the main differences between 
the configurations, as one can deduce by comparing the 
left (native-state) and right (activated-state) part with 
respect to the diagonal in Figs, la and lb. Here, we 
have found that the most significative contrast between 
the native and activated state is obtained for R =6 A in 
bR and for R =12 A in rat 0R-I7. 

The graph representation is then used to describe the 
protein physical features. The graph turns into the 
INPA, by attributing an elemental impedance Zij = 
h,j/Aij{p~^+iei^jeoUj) [lol[Tl| . to each link. The elemen- 
tal impedance is constituted by a an Ohmic resistance 
connected in parallel with a parallel plate capacitor, filled 
with a dielectric, so that both charge motion and elec- 



trical polarization can be described. Finally, the model 
was implemented to determine non-linear I-V character- 
istics by adopting a sequential tunneling mechanism of 
charge transfer, as detailed in Refs. (|l3h . We assume 
a tunneling probability as given by the WKB approxi- 
mation: Pi j = exp[— ^^•\/2m($ — eVi,j)], where Vij is 
the potential drop between the i-th and j-th node, m 
is an effective electron mass, here taken as that of the 
free electron, and $ is the barrier height, taken, for these 
calculations, equal to 59 meV. 

Figure 2a reports the Nyquist plots for the native and 
activated state of bR and for R =6 and 12 A , respec- 
tively. The plots exhibit an almost semicircular shape, 
slightly suppressed because of the presence of many 
relaxation-time scales [l^ . The contrast between the two 
states takes the highest level with R =6 A , lowering sys- 
tematically for increasing cut-off values. In particular, 
the static impedance is found to decrease for about 10%, 
at most, when moving from the native to the activated 
state. To complete the analysis of the impedance spec- 
trum, figure 2b reports the frequency dependence of the 
real, Re{Z {uj)} , and imaginary, Im{Z(u))}^ parts of the 
impedance (Bode plots), respectively, for the highest con- 
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FIG. 4: (Color online) I-V characteristics of the native and 
activated states of rat 0R-I7. Empty circles refer to the native 
state, and full circles to the activated state, respectively. The 
continuous lines are guides for eyes. 



trast of R =6 A . To our knowledge, experimental data 
on bR are not available. 

Figure 3a reports the Nyquist plot for the native and 
activated state of rat 0R-I7 and for R =6 and 12 A , 
respectively. In analogy with the case of bR, the plots 
exhibit an almost semicircular shape. The contrast be- 
tween the two states takes the highest level with i?=12 
A . In particular, the static impedance is found to de- 
crease for about 60%, at most, when moving from the 
native to the activated state. Both the shape of the plot 
and the magnitude of the resistance charige are in qual- 
itative agreement with EIS experiments [6| obtained for 
a self assembled layer of 0R-I7 in the presence of its spe- 
cific odorants heptanal and octanal. In Fig. 3a, the two 
points under the arrows on the x axes report the dose 
response of 0R-I7 layered structure to a concentration 
of 10"''^ mol/L of octanal and heptanal. In analogy with 
the case of bR, Fig. 3b reports the normalized real and 
imaginary parts of Z (Bode plots) as function of w, for 
the highest contrast of i? = 12 A , respectively. We would 
remark the difference between the Nyquist plots of bR 
and 0R-I7, the former showing the maximal resolution 
of the native and activated configurations for i? = 6 A , 
the latter for i? =12 A . The origin of this difference is at- 
tributed to a more massive conformational change in the 
olfactory receptor than in the bR, as already documented 
by the contact maps in Fig. la and lb. 

Figure 4 reports the I-V characteristics for the rat OR- 
17, obtained with R =6 A and the same transport param- 
eters used for bR [ll] . With respect to bR, we have found 
a larger difference between the current response of the na- 
tive and activated state, with the native state showing a 
quasi-linear response in the range of considered voltages. 
By contrast, the current corresponding to the activated 
state exhibits a significant super-linear increase, in close 



analogy with what observed in bR experiments. This 
result can be interpreted in terms of large distances be- 
tween amino-acids in the native state, which oppose the 
tunnelling mechanism of charge transfer, and of a sig- 
nificant decreasing of the same distances after the con- 
formational change, which favors charge transfer. This 
interpretation applies well also to the large difference ex- 
hibited by the shapes of the Nyquist plots for the 0R-I7, 
and confirms the strict correlation between the protein 
tertiary structure and its conformational change. 

Finally, we stress that present theoretical results take 
into account only the topological structure and its mod- 
ifications, and refer to a protein completely activated. 
Furthermore, different bio-chemical contributions, as the 
role of the environment, are neglected. Indeed, be- 
ing independent of configurational change, the environ- 
ment should play a minor role in determining the change 
of electrical properties induced by the conformational 
change, which justifies the omission of its effects. 

In conclusion, we have investigated the correlation be- 
tween a conformational change and the variation of elec- 
trical transport properties in a sensing protein. Theory 
makes use of an appropriate procedure devised to account 
for the interaction responsible for charge transfer between 
neighbouring amino-acids. On this basis, we have inves- 
tigated the change of the small-signal impedance spec- 
trum and that of the static I-V characteristics in bacteri- 
orhodopsin and rat 0R-I7. Concerning the small-signal 
impedance spectrum, the theory is compared with EIS 
experiments carried out on self-assembled monolayers of 
0R-I7, an odour sensitive GPCR protein, over function- 
alized gold substrates. Qualitative agreement with exper- 
iments is found when considering a maximum interaction 
radius between amino acids of 12 A , which evidences a 
substantial modification of the Nyquist plot (decrease of 
the resistance up to a maximum of about 60%) in the 
presence of specific odorants, like heptanal and octanal. 
Concerning the I-V characteristics, the theory applied to 
a nanolayer of 0R-I7 evidences a substantial increase of 
the current, at a given voltage, when going from the na- 
tive to the activated state. By testing the predictability 
of the model, for the Nyquist plot of bR we found a small 
but still significant deviation of the spectrum with a de- 
crease of the resistance for a maximum of about 10 % 
when going from the native to the activated state. 
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